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ABSTRACT 

We present evidence for widespread relativistic effects in the central regions of 
active galactic nuclei. In a sample of 18 Seyfert 1 galaxies observed by ASCA, 14 
show an iron Ka line which is is resolved, with mean width a-^ a = 0.43 ± 0.12 keV 
for a gaussian profile (Full Width at Half Maximum, FWHM~ 50,000 km s" 1 ). 
However, many of the line profiles are asymmetric. A strong red wing is indicative of 
gravitational redshifts close to a central black hole and accretion disk models provide 
an excellent description of the data. 

The peak energy of the line is 6.4 keV, indicating that it arises by fluorescence in 
near-neutral material. Our fits imply a low inclination for the disk in these Seyfert 1 
galaxies, with a mean of 30°, consistent with orientation-dependent unification schemes. 
Differences in the line profiles from source-to-source imply slight variations in geometry, 
which cannot be accounted for solely by inclination. In most cases, we require that 
the line emission arises from a range of radii. Although a small contribution to the 
emission from a region other than the disk is not ruled out, it is not generally required 
and has little effect on our conclusions regarding the disk line. Our data are fit equally 
well with rotating (Kerr) and non-rotating (Schwarzschild) black hole models. We 
find a mean spectral index in the 3-10 keV range of < T3_io >= 1.91 ± 0.07 after 
accounting for the effects of reflection. 

Such observations probe the innermost regions of AGN, and arguably provide the 
best evidence yet obtained for the existence of super-massive black holes in the centers 
of active galaxies. 

Subject headings: galaxies: active - galaxies: nuclei - galaxies:Seyfert - X-rays:galaxies 

1. Introduction 

The X-ray spectra of Seyfert 1 galaxies can contain numerous components. The continuum 
has an apparent power law form, which covers several decades in energy. This is attenuated at 
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soft X-ray energies by absorption both in the interstellar medium in our Galaxy and by material 
local to the active galaxy. An emission line from iron Ka has been observed at 6-7 keV in the vast 
majority of sources. A hard component affects the spectrum above ~ 10 keV which is probably 
produced by Compton scattering ("reflection") of the continuum in the material which produces 
the line. In the soft X-ray regime, many sources have "soft excess" emission and can show emission 
features from oxygen and/or iron-L. 

At least three nuclear emission/absorption regions are envisaged to account for these 
components. The region which emits the primary X-ray continuum may consist of hot and/or 
relativistic electrons and/or electron-positron pairs which Compton up-scatter softer "seed" 
photons. Compton-thick material in a low state of ionization, which may be identified with 
an accretion disk, gives rise to the iron Ka line, the reflection continuum and possibly the 
"soft excess". In many cases, the local absorption is due to a Compton-thin "warm absorber": 
highly-ionized gas identified by features due to oxygen and iron, which may also produce the soft 
X-ray lines. 

Detailed measurement of the various emission and absorption features, and continua, can 
therefore yield information regarding the physical conditions of the regions, such their dynamics, 
geometry and ionization state. Of particular interest in this regard is the iron Ka emission 
line. These lines were first discovered with HEAO-1 in heavily obscured Seyfert galaxies, such 
as NGC 4151 (Warwick et al. 1989 and references therein), and a number of Narrow Emission 
Line Galaxies (Mushotzky 1982). Not until the launch of Ginga was iron Ka emission recognized 
as an important, and possibly universal, property of Seyfert galaxies (Nandra & Pounds 1994, 
hereafter NP94, and references therein). From the time of their discovery in relatively unabsorbed 
AGN, the emission lines have been assigned an origin in material intimately associated with the 
accretion process, the most promising candidates being an accretion disk, or "blobs" of material 
surrounding the central source (Guilbert & Rees 1988; Fabian et al. 1989; Nandra et al. 1989; 
Pounds et al. 1990). However, for individual sources, the properties of the line were not generally 
well determined, given the limited spectral resolution of Ginga. 

NP94 have discussed the properties of the iron Ka line in some detail, using a sample of 
Seyfert galaxies. The mean energy was found to be ~ 6.4 keV, which indicates that the most 
likely origin for the line is fluorescence in near-neutral material. Given this information, it can be 
deduced that the line cannot arise in material uniformly covering the source. The mean equivalent 
width of the line of 140 eV requires a high optical depth for the material, which would produce 
considerably more soft X-ray absorption in these objects than is observed. Indeed, the equivalent 
width of the line is sufficiently high to indicate that the material may be optically thick to 
Compton scattering. Further evidence for a Compton-thick region exists in the form of the "hard 
tail" , which can be produced by "reflection" from the line-producing region. Given this evidence, 
an accretion disk is an obvious source for these features, especially when considering that the 
strengths of both the line and "reflection" component imply a covering fraction of ~ 50 per cent 
for the material. Nonetheless, the "blob" geometry was not ruled out by the Ginga data (Bond & 
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Matsuoka 1993; Nandra & George 1994). 

Further impetus for the study of iron Ka lines in Seyfert galaxies has come from more recent 
ASCA data, which have indicated that this feature may be a key element in our understanding of 
AGN. The earliest data showed good evidence that the emission lines were resolved (Fabian et 
al. 1994; Mushotzky et al. 1995), confirming tentative suggestions from the Ginga data (NP94). 
Most dramatically, high signal-to-noise data for two sources have shown characteristic line profiles 
(MCG-6-30-15 Tanaka et al. 1995; NGC 4151 Yaqoob et al. 1995). The lines in these sources are 
extremely broad, with FWHM implying relativistic velocities of order 0.2c. Furthermore, there is 
a strong asymmetry to the red, which is indicative of the gravitational redshifts associated with 
the inner regions of an accretion disk surrounding the black hole. Fabian et al. (1995) concluded 
that this was the most plausible explanation for such profiles, rejecting other interpretations such 
as Comptonization. Whilst it is tempting to extrapolate such a model and apply it to all AGN, we 
must bear in mind that these extraordinary profiles are relatively rare in the literature. Further 
reinforcement of the black hole/accretion disk paradigm requires confirmation of that model for 
further sources. It is our intention to explore such avenues in the present paper. 

In a previous paper (Nandra et al. 1996b; hereafter paper I), we have presented imaging and 
timing data for a sample of Seyfert 1 galaxies observed by ASCA. Here we present corresponding 
spectral data in the 3-10 keV regime in order to define the properties of the iron line. In Section || 
we describe our analysis techniques. We will present the spectra in the full band in a subsequent 
paper (hereafter paper III). Following that, Section |3] present simple parameterizations for the 
iron Ka line. In Section |4] we test specific disk-line models. Our aim thereby is to determine if 
such models, which have successfully described the two sources mentioned above, are generally 
applicable. In Section ||, we compare our results to previous X-ray observations. We summarize 
our findings in Section ^ and discuss them in Section [?]. 



2. Data analysis 

In the present paper, we examine the spectra of a sample of Seyfert 1 galaxies observed by 
ASCA in the 3-10 keV band. We chose this energy range as it is relatively free from the effects of 
absorption, either by cold gas or the ionized "warm absorber", allowing us to attempt a definition 
of the continuum shape and, particularly, the iron Ka emission line. Furthermore, we consider 
only the time-averaged spectrum for each observation, to maximize the signal-to-noise ratio. Our 
sample is listed in Table We refer to each observation using the name of the source and, in cases 
where there are multiple observations, the observation number in brackets (e.g. NGC 4151(2) is 
the second observation of NGC 4151; note that rejected observations are still assigned a number 
- see paper 1 for details). The selection criteria, observation log and data analysis methods are 
presented in Paper 1. Briefly, we have employed data from the US public archive at NASA/GSFC. 
We applied standard criteria for rejecting poor quality data and then extracted the source spectra 
using the FTOOLS/XSELECT package from both the Solid-state Imaging Spectrometers (SIS) 
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and Gas Imaging Spectrometers (GIS). Background spectra were extracted from source-free 
regions. 

For both instruments we considered only the gain-corrected (Pulse Invariant or PI) channels. 
The raw spectra were rebinned such that each resultant channel had at least 20 counts per bin, 
permitting us to use x 2 minimization for spectral fitting. All such fits were undertaken using the 
XSPEC package (v9.0). Calibration data were taken from the HEASARC calibration database. 
The SIS response matrices were released on 1994 November 9 and the GIS matrices on 1995 
March 6. Effective areas were calculated using the standard software, which does not include a 
parameterization of the azimuthal dependence of the the point spread function. This can cause 
inaccuracies in the normalization, particularly in the case where the source region is not circular. 
Uncertainties in determining the source centroid can also result in normalization errors, and the 
cross-calibration of the four instruments is not known to arbitrarily high accuracy. These factors 
lead to cross-calibration uncertainties between the four instruments, but as the response of the 
XRT is not strongly dependent on energy, we do not expect them to introduce distortions into 
the individual spectra. We therefore chose to fit all four ASCA detectors simultaneously with the 
same model, but allowed the normalization of that model to be free for each detector. 

Due to radiation damage, the response of the CCD detectors is changing slowly with time 
(Dotani et al. 1995). The main effects of this damage relevant to the analysis presented here 
are small changes in the energy scale and a degradation of the energy resolution. We have not 
accounted for such effects, which are most pronounced in 4-CCD mode data. However, even in 
this most extreme case, we expect the effects to be small for our data, which were taken early 
in the mission. The estimated offset and resolution changes are < 50 eV for 4-CCD mode and 
are negligible for 1 and 2-CCD modes for data taken < 1 yr after launch (Dotani et al. 1995). 
Normalizations and fluxes quoted for the sources will be for the SISO detector. Typically, the two 
SIS detectors agree in flux to within ~ 1 — 3 per cent, but are higher than the GIS by ~ 10 — 20 per 
cent. 

2.1. Statistical considerations 

As mentioned above, we have employed x 2_mm i m i za tion for spectral fitting (Lampton, 
Margon & Bowyer 1976). Although this forces us to rebin the data arbitrarily, sometimes beyond 
the spectral resolution of the instrument, we prefer this over maximum likelihood methods such 
as the C-statistic (Cash 1979) because the data can be background-subtracted. This is highly 
desirable in cases, such as these, where reliable background models have not yet been developed. 
Although the background represents typically only a few per cent of the total count rate in the 
3-10 keV band, it can have a significant effect at high energies. We therefore note here the fact 
that we have not accounted for vignetting of the diffuse X-ray background when calculating our 
background rates, nor any other variations in the background rates over the detectors. Such 
effects are not currently well-understood for ASCA data. Our chosen bin-size is not optimal 
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for determining the goodness-of-fit; at low energies the data are under binned compared to the 
resolution of the detector. We prefer this to the alternative of over binning, which would sacrifice 
spectral information essential to our determination of the line parameters. Although detailed 
justification for our choice of 20 counts per bin requires simulations beyond the scope of this paper, 
a comparison of our spectra with the simulations of Nousek (1990), suggests that using a gaussian 
approximation and x 2 - m i n i m i z ation will introduce uncertainties of the order of 1 per cent to our 
derived parameters, smaller than the typical statistical error. 

Once the minimum x 2 value is found in each fit, confidence regions for the parameters have 
been determined using the prescription of Lampton et al. (1976). Unless otherwise stated we have 
calculated the confidence region by projection of the joint confidence region defined by the number 
of "interesting" parameters, as recommended by Lampton et al.We have considered the model 
normalizations as uninteresting. This is because we do not expect the origin of the difference in 
the normalizations - uncertainties in the effective areas - to affect the spectral form. Thus the 
number of interesting parameters which we have used to calculate the \ 2 deviates is given by 
Nf r — 4, where iVf r is the number of free parameters in the model. Specific A% 2 values are shown 
as footnotes to the Tables. For completeness we have quoted the best-fit normalization for one of 
the models in Table although we have not attempted to determine the confidence interval for 
that parameter, as we consider it uninteresting. 

In most cases we quote the 68 per cent confidence limits for the parameters, which is 
the same confidence level as the lcr distribution of a gaussian. Furthermore, in some specific 
circumstances, we have gone on to interpret the 68 per cent confidence limit as if it were the la 
error bar of a gaussian distribution. Examples are the calculation of the weighted mean, and the 
determination of the mean and dispersion using the maximum-likelihood method. This assumes 
that the confidence surface is itself gaussian, which is often a poor approximation. However, this 
is unimportant, as long as the gaussian approximation does not underestimate the true confidence 
region. We believe this to be true in most circumstances. However, in most cases the confidence 
region is asymmetric, resulting in two-sided "error bars". In all such cases we use the larger 
uncertainty as our la error bar. Furthermore, where the 68 per cent confidence region overlaps the 
limits of the allowed parameter values (e.g., where the inclination reaches 0°or 90°), we assume 
that the la error bar is essentially infinite, thus placing no weight on the derived value. This is 
a highly conservative approach which should avoid any spurious results arising from our use of 
confidence regions as gaussian error-bars. 

In comparing various models for the X-ray spectra, and particularly the iron lines, we have 
employed various statistical tests, including: the acceptability of x 2 > although as noted above, 
our data are not optimally binned to maximize x 2 i the F-test (e.g., Bevington 1969) which has 
been employed to determine whether additional components, or parameters, provide a significant 
improvement; the likelihood ratio (Edwards 1972; Mushotzky 1982) which has been used to 
determine whether models with the same number of free parameters are more or less likely to 
produce the observed data. We have adopted the 95 per cent confidence level in assessing whether 
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or not various features are real, but we bear in mind that for a sample of 23 (or more often 22) 
observations, this will typically result in one spurious detection each time the sample is considered 
as a whole. 

3. Iron Ka emission: simple parameterizations 

We first attempted some simple spectral fits to the ASCA data, without reference to a specific, 
physical model. 

3.1. Narrow lines 

Our initial fit was of a power-law with a gaussian component with width, a = 10 eV, to 
represent the iron line. This width is smaller than the instrumental response and therefore is 
effectively "narrow" for our purposes. In practice, this fit will model any line emission up to the 
width of the spectral response of ~ 120 eV (FWHM) at 6 keV. The gaussian provides a significant 
improvement in the fit in all 18 sources, confirming that iron Ka emission is extremely common in 
Seyfert 1 galaxies. The details of these narrow-line fits are shown in Table The energies quoted 
are in the rest-frame throughout. Three of the fits are unacceptable at > 95 per cent confidence: 
NGC 4151(2), NGC 4151(5), MCG-6-30-15(l). The mean values of the parameters derived from 
these fits are shown in Table ||. That table shows the unweighted mean and standard deviation 
(column 1), the weighted mean and uncertainty (column 2), and in columns 3 and 4 estimates 
of the mean and dispersion of the parent population, calculated using the maximum likelihood 
method of Maccacaro et al. (1988). We consider the latter two values to be most appropriate 
for determining the average properties of the sample, and quote these in the text. A histogram 
of the line energies and equivalent widths is shown in the upper panels of Fig. [I]. Confirming 
and extending the Ginga results, our ASCA observations show a strong preference for iron in a 
relatively low state of ionization, with the mean energy of < E-& a >= 6.36 ± 0.02 keV, which 
is even indicative of a small net redshift (but only at ~ 90 per cent confidence). There is no 
measured dispersion in these values, with an upper limit of 60 eV; the large deviation in the case of 
NGC 6814 is due its large measurement error, as it is a factor > 10 weaker than all other sources. 
The lack of dispersion is remarkable; any systematic deviations in the energy scale due to residual 
calibration errors (see Section |2|) must be smaller than this. Also note that the differences in the 
observed- frame energy between the sources is ~ 300 eV. We find a mean ionization state <Fe xvi 
(-Exa < 6.41) at the 99 per cent confidence level. The mean equivalent width of 98 ± 12 eV for 
these narrow lines is only ~ 70 per cent of the corresponding value derived by Ginga (140 eV; 
NP94). An examination of the residuals makes clear why this is the case. In a number of sources 
the line flux clearly extends beyond the bounds of the spectral response function, indicating 
the line has a significant width. The wider response of Ginga would result in more of this flux 
being modeled in a fit what what is effectively a delta-function profile. Such a profile is clearly 
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not appropriate for our higher-resolution ASCA data, and we now go on to describe this more 
quantitatively. 

3.2. Broad lines 

Some Ginga spectra of Seyfert galaxies showed a significant improvement in the fit when the 
line width was allowed to be a free parameter (NP94). Due to the nature of the Ginga spectral 
response this implied huge line widths of order 1 keV in the cases where significant broadening was 
indicated. This led NP94 to speculate that the apparent broadening might be spurious, instead 
reflecting unmodelled continuum components. Improved spectral resolution is the key to removing 
this ambiguity. Indeed, early ASCA observations showed that, at least in some cases, the iron 
line was significantly broadened (e.g. Mushotzky et al. 1995). We tested this for our sources by 
allowing the width of the gaussian a to be a free parameter. The results of these fits are shown in 
Table g 18/23 observations and 14/18 sources showed a significant improvement, with only two 
fits just unacceptable at the 95 per cent confidence level: NGC 4151(2) and MCG-6-30-15(l). In 
cases where the line width is not detected, we have shown only the 90 per cent upper limit to the 
line width. These are all consistent with the mean value of < ctkq >= 0.43 =L 0.12 keV (Table [|). 
This corresponds to a FWHM of approximately 47,000 km s _1 . A histogram of the line widths is 
shown in Fig. 

Once again, the energy of these broad lines strongly indicates material in a low ionization 
state is the source of the line, with a mean energy < E-^ a >= 6.34 ± 0.04 keV. The 99 per cent 
confidence upper limit to the energy of 6.41 keV again corresponds to an ionization state for 
iron <Fe xvi - above this ionization state the line rises sharply in energy (Makishima 1985). 
If the material were in collisional ionization equilibrium, this corresponds to a temperature of 
<^ 3 x 10 6 K (Arnaud & Raymond 1992). Of course, any material close to the central source is 
more likely to be photoionized than collisionally ionized, in which case the temperature is probably 
much lower. 

The mean equivalent width of these broad lines is < W~K a >= 160 ± 30 eV. This value is 
consistent with the mean determined by Ginga for a line of a = 0.1 keV, of 140 ± 20 eV. As 
mentioned above, a fit with a line which is narrow compared to the instrumental response of a 
given instrument will tend to model any line flux which has a width up to the spectral resolution. 
For the case of Ginga this is equivalent to a ~ 0.5 keV (Turner et al. 1989). This corresponds 
roughly to the mean width of the lines determined in our broad-gaussian fits. In this respect, then, 
the equivalent widths determined by ASCA and Ginga show a remarkable consistency. 
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3.3. Line Profiles 

At least in some sources, the iron Ka line is not well modeled by a symmetrical profile such 
as a gaussian. We have attempted to determine the line profiles in a relatively model-independent 
way, with the results show in Fig ||. These were determined with a fit to the 3-10 keV SIS data only 
with a power law, but excluding the iron band between 5-7 keV. We then multiplied the resultant 
data/model ratio by the best-fit continuum model. This results in an estimate of the line profile 
as seen in the observed frame, deconvolved from the instrumental response. Almost all of the 
sources show a clear excess, with broadening far greater than the instrumental response. It is also 
evident that a number of sources show skewed, asymmetric profiles. This could represent multiple 
line components (for instance, there may be both broad and narrow components). However, many 
sources show the characteristic black hole/accretion disk profile (Fabian et al. 1989; Laor 1991; 
Matt et al. 1992), in which case they may be predominantly one component. 

To demonstrate this further we have created a composite line profile for all 18 sources. This 
was achieved from by transforming the above-mentioned data/model ratios into the rest frame of 
each source and taking the mean ratio in bins of width 40 eV. To estimate the mean line profile 
from these plots, we have multiplied by a continuum model defined by the mean values for the 
sample. This is shown in Fig ||a. In Fig ^|b we show the same plot excluding MCG-6-30-15 and 
NGC 4151, to investigate whether these sources dominate the mean profile. Several things are 
apparent from Fig. 4: 

• The line profiles are clearly consistent, showing that our results are not biased by the 
inclusion of MCG-6-30-15 and NGC 4151. 

• The profiles are extremely broad ~ 2 keV Full-width at Zero Intensity (FWZI), which 
corresponds to velocities of order 0.3c. 

• As in many of the individual sources the profile consists of a relatively narrow core, with an 
underlying broad component. 

• The core peaks at an energy remarkably close to 6.4 keV. 

• The broadening is primarily to the red, with relatively little flux blue-ward of 6.4 keV and 
specifically no strong component due to helium-like or hydrogen-like iron (6.7-6.9 keV). 

To illustrate these more quantitatively, we have modeled the profiles with two gaussians, 
following Tanaka et al. (1995). The resulting fit is shown in Fig. |I]a. The centroid energies of the 
core are 6.38 keV in both cases, with a width of a = 0.1 keV. The broad component is strongly 
redshifted with centroid energy ~ 6 keV and width of a = 0.7 keV and carries ~ 75 per cent of 
the total flux. The dotted line in Fig. |I]b shows the instrumental response. Whilst some artificial 
broadening could result from radiation damage to the CCDs, particularly in 4-CCD mode, these 
residual calibration uncertainties of ~few per cent in this regime are much smaller than the 
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observed deviations (~ 10 — 20 per cent in the red wing). To demonstrate this, we show in Fig. 
the residuals for our Seyfert sample compared to those for the Coma Cluster and also the Dark 
Earth (essentially the particle background of the SIS; Gendreau 1995). All three plots show an 
emission line, but only the Seyferts show the red wing, indicating that this feature is not due to a 
persistent calibration error. 

It is also highly unlikely that emission lines from any any other element contribute significantly 
to the emission in the 5-7 keV band. The other elements with K-shell energies (V, Cr, Mn, Co) in 
this range have abundances and fluorescence yields which would lead to line intensities of a factor 
> 100 less than iron, to which we attribute the entire line flux. Iron line blends alone cannot 
account for the observed widths in most cases, as the line flux extends up to ~ 1.5 keV below the 
energy expected from neutral iron. Compton down-scattering can, in principle, redshift the line 
sufficiently, but many workers have argued against this mechanism (Mushotzky et al. 1995; Tanaka 
et al. 1995; Fabian et al. 1995). The dramatic redshifts and broadening expected due to special 
and general relativistic effects close to a black hole provide a plausible explanation for the observed 
line profiles. One class of models which has been successfully applied to ASCA observation of 
MCG-6-30-15 and NGC 4151 is that of an accretion disk (Tanaka et al. 1995; Yaqoob et al. 1995). 
The other sources apparently exhibit rather similar profiles and we now consider such models. 

4. Accretion Disk modeling 

We will refer to a number of parameters in the following sections which characterize the disk 
line models: the inner radius of the disk, R\, the outer radius, R Q , the inclination of the axis of 
rotation with respect to the line-of-sight, i, the rest energy i?j( a and the line normalization /kq • 
We excluded NGC 6814 from the following analysis as the signal-to-noise ratio is so much lower 
than that of the other sources, leaving 22 observations of 17 sources. 

4.1. Schwarzschild model 

Our first test of the black hole/accretion disk hypothesis was using the model of Fabian et al. 
(1989). Their disk-line model assumes a Schwarzschild geometry and computes the line profile 
only and not the strength. No specific geometry for the X-ray source is assumed, and therefore the 
line emissivity is parameterized by a power law as a function of radius, R~ q . The value of q can 
be a free parameter. For a point-like X-ray source located above the center of the disk, we expect 
q ~ in the inner regions, steepening to q ~ 3 at large radii. However, different emissivity profiles 
may be relevant, for example if there are multiple X-ray sources, or if the source is not point-like 
compared to the line-emitting region. 

The line profile has a strong dependence on q, as for values > 2, the emission is concentrated 
in the inner disk, where extreme gravitational and Doppler effects operate. In these circumstances, 
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R Q can rarely be constrained. Conversely, if the outer disk is the dominant source of line emission, 
poor or no constraints can be placed on R^. We have therefore chosen to fix the inner and outer 
radii of the disk at 6 R g and 1000 R g respectively (where R g is the gravitational radius of the 
black hole) and have used q as a parameterization not only of the geometry of the X-ray source, 
but of the accretion disk itself. We allow the inclination of the disk to be a free parameter and also 
the normalization of the line. We will compare the strength of the lines to the disk model later. 
Finally, we assume Eka = 6.4 keV in the rest frame. As we showed in Tables || and ||, and Figs || 
and || both the peak and mean line energies lie close to this value in all cases, with no evidence for 
a helium-like or hydrogen-like line. 

The results of these disk-line fits are shown in Table There is a substantial improvement in 
the fit compared to the broad gaussian model in a large number of cases, due to the asymmetry of 
the profiles. All of the disk line fits are acceptable at < 95 per cent confidence. The total x 2 f° r an 
22 observations was 16762.4 for the broad line fit, and 16627.3 for the disk lines (Ax 2 =135.1) which 
is highly significant as measured by the likelihood ratio (see Section Q), which is shown in Table |8|. 
The mean power-law index is similar to the broad-line fits < T3_io >= 1.78 ± 0.07. However, 
the line equivalent width is rather larger, at < Wk q >= 290 ± 50 eV. Typically, a low inclination 
(~face-on) disk is preferred in our fits. The mean inclination angle is < i >= 29 ± 3 degrees. If 
the accretion disk were oriented randomly, we might expect a mean inclination of ~ 60°, which is 
not preferred in our fits. We comment on this later. The mean value of q, which parameterizes the 
geometry of the system is < q >= 2.5 ± 0.4. This is in good agreement with the value expected 
for a point-like X-ray source above a flat disk, as discussed by Matt, Perola & Piro (1991), and is 
therefore supportive of the disk model. Histograms of i and q are shown in Fig. ^. 

Likelihood analysis reveals no evidence for a significant dispersion in the values of q obtained 
in our fits, which suggests a common geometry for the sources. However, our conservative error 
prescription, and the assumption that the 68 per cent confidence limit equates to a la error bar 
makes it difficult to rule out such an hypothesis based on this analysis. A more sensitive method of 
determining whether or not there are differences in geometry from source-to-source (in addition to 
differences in inclination) is to fit the data explicitly with q fixed at the mean value, and compare 
the x 2 value so obtained with that when q was left free. We have performed such fits, and these 
are tabulated in Table || It is clear from column 5 of Table [|, which shows the F-statistic for the 
addition of q as a free parameter, that many sources prefer a value of q different to the mean. 
Indeed 11 of 22 observations and 8/17 sources show an improvement at > 95 per cent confidence. 
As an ensemble, the reduction in x 2 is 96.1 for the addition of 22 additional free parameters, 
significant at > 95 per cent confidence. This implies that for this model, the differences in line 
profiles cannot be accounted for solely by difference in inclination; assuming that the whole of the 
emission line arises in an accretion disk, we infer that the X-ray source/accretion disk geometry is 
different from source to source. 

The above model assumes that the line profile consists of emission integrated over a range of 
radii, each of which gives rise to a characteristic line shape. Variations in q correspond to different 
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weightings over the disk. An alternative explanation for the range of observed profiles is that 
the emission is dominated by a single radius in each source, but that this characteristic radius is 
different in each case. In order to investigate this possibility we have repeated the disk line model, 
but setting R\ = R Q and allowing that radius to be free, q is clearly irrelevant in this case. The x 2 
values for these fits are shown in Table |8[ For the sample as a whole they are considerably worse 
than those with a range of radii with an increase in x 2 °f 106.5 for the same number of degrees of 
freedom. In only two observations is there an appreciable improvement (i.e., a significant reduction 
in x 2 ), NGC 4051 and MCG-6-30-15(l). 

4- 1.1. Rest energy of the line 

In the presence of substantial gravitational and Doppler effects inferred above, the energy 
derived from a gaussian fit may not reflect the true rest energy of the line. For example, it is 
conceivable that gravitational redshift might cause a line which was originally helium-like, with an 
energy of 6.7 keV, to appear at a lower energy. Given the remarkable consistency of our gaussian 
fits with an energy of 6.4 keV and the lack of any strong emission above that energy (Fig. this 
seems rather contrived. As shown by Matt et al. (1992), for an ensemble of disk lines, the mean 
energy always lies close to the rest energy, for a variety of emissivity laws. This is due to the fact 
that the peak energy is dominated by the low- inclination profiles, which are of highest intensity 
(George & Fabian 1991). We have investigated this further by constructing summed model profiles 
for an ensemble of accretion disks with inclinations distributed uniformly in cosi. For q < 3, the 
peak energy for the ensemble is redshifted by < 1 per cent of the rest energy. Therefore, unless the 
emission originates entirely within the central regions, the peak energy in Fig || strongly indicates 
a near-neutral origin for the emission line. This is further illustrated in Fig. which shows the 
summed model profile for the disk line fits shown in Table |j. 

Nonetheless, we have gone on to explore whether our observed profiles are consistent with a 
helium-like line distorted by the gravitational and Doppler effects in the inner disk. Our first test 
was with a model identical to that used in the fits of Table |3], i.e. R\ = 6, R Q = 1000, q and i free, 
but with the energy of the disk line fixed at 6.7 keV in the rest frame. These fits were considerably 
worse than with a 6.4 keV line. The total x 2 was 16737.2 for 16973 d.o.f., compared to 16627.3 
in the near-neutral case. We infer the latter to be > 1000 times more likely (Table ||). We have 
investigated whether this might be due to our assumption of R Q = 1000. At large radii, the lack 
of any significant gravitational or Doppler shifts results in emission close to the rest energy of the 
line, yet we observe little emission around 6.7 keV (Fig. § Fig. ||). Relaxing the constraint on R Q 
with the helium-like disk line does improve the fits. However, these fits were still slightly worse 
(total x 2= 16633.2), than those with the 6.4 keV line and a fixed outer radius (total x 2= 16628.2), 
despite the additional 22 free parameters. However, examination on a case-by-case basis suggests 
that 3C 120 shows the most marked difference in x 2 - Removal of that outlier would leave the 6.7 
keV/free R Q fits showing a slight improvement over those at 6.4 keV/fixed R Q , but not a significant 
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one. 

As a further test, we have allowed R Q to be free in the 6.4 keV fits, for comparison with the 
helium-like model. Our suspicion that R D could not be constrained by these data was largely 
confirmed by this analysis, with only four observations, NGC 4051, NGC 4151(2), MCG-6-30-15(l) 
and Mrk 841(1), showing a significant improvement at > 95 per cent confidence. The lack of other 
detections is likely to be due primarily to signal-to-noise ratio, the degeneracies inherent in the 
disk line models and the fact that the bulk of the line emission is produced in the central regions. 
The total x 2 f° r those fits was 16597.0, indicating that a 6.4 keV line is > 1000 times more likely 
(Table ||) than a 6.7 keV line, even if the outer radius is left free. That conclusion holds even if 
3C 120 is removed from consideration. 

Relaxing the constraint on the inner radius is unlikely to improve the fits which assume 
helium-like iron, as we have fixed the inner radius at the last stable orbit in the Schwarzschild 
metric. This maximizes the gravitational effects, which would tend to bring any highly ionized 
emission closer to the observed values. 

For the four line models tested thus far (narrow gaussian, broad gaussian, disk line with fixed 
R Q and disk line with free i? Q ), we have shown that a 6.7 keV line is strongly disfavored. However, 
this analysis does not preclude the possibility that more complex models can be constructed which 
conspire to produce an apparent rest energy of 6.4 keV. Particularly, we have assumed that the 
line can be described by a single component. We now explore the possibility that an additional, 
narrow line is present. 

4-1.2. Narrow line contribution 

It is possible that there is a contribution to the emission line from another region. Differing 
relative contributions from that region and the disk could then confuse our results from the disk 
line fits. A promising candidate is the obscuring torus hypothesized in Seyfert 1/2 unification 
schemes, which could produce a significant contribution in some circumstances (Krolik & Kallman 
1987; Ghisselini, Haardt & Matt 1993; Krolik, Madau & Zycki 1993). A relatively narrow line 
might also be produced from the optical "broad line" region or even the warm absorber. 

We have tested for the presence of such a feature by adding a narrow gaussian, at a rest 
energy of 6.4 keV to the disk-line models described in Table This provided a significant 
improvement in the fit for 6 observations: NGC 3516, NGC 3783(1), NGC 4151(2,4,5) and NGC 
5548. The improvement for the ensemble is Ax 2 =68.9 for the 22 degrees of freedom, which is 
formally highly significant. The improvement is primarily due to NGC 4151 - without that source 
the improvement for the sample would not be significant (A% 2 =32.3 for 19 degrees of freedom). 
This is noteworthy, as NGC 4151 is the only heavily-absorbed source in our sample; conceivably 
the line-of-sight column could be contributing to the emission line in this case (see Weaver et al. 
1994 and Yaqoob et al. 1995 for detailed discussion of the emission line in NGC 4151). In the 
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other cases where the narrow line provides a significant improvement, this may simply reflect the 
fact that we have over-restricted the parameter space for the disk-line model. Apparently there is 
no requirement for a narrow line contribution in the majority of our sources, although the upper 
limits are not restrictive (typically ~ 100 eV, consistent with predictions for the torus; Ghisselini 
et al. 1993). The mean equivalent width of the narrow component for the sample is only ~ 30 eV. 

Even if there is such a narrow component, it has little effect on our conclusions regarding the 
remainder of the line. This is illustrated in Fig. which shows the summed model profiles with 
and without the narrow line. To quantify the changes we show in Fig. || the change in the disk 
line parameters (Wk , i, q) when the narrow line is added to the model. Obviously, the equivalent 
width of the disk line decreases when the narrow line is added, but typically only by a few tens of 
eV. There effect on the inclinations is negligible; our fits are driven to low inclinations by the lack 
of significant blue flux, which is not changed by the narrow component. There is a small effect 
on the emissivity index, with the the mean value of q from the disk-line-plus-narrow-line fits was 
q = 2.8, slightly steeper than that without the narrow line, but not significantly so. We have also 
explored whether the presence of a narrow line can explain the apparent difference in geometries 
implied by the range of q values in Table ||. 

Fixing q at the mean value of 2.8 resulted a significantly worse fit in 9/22 of the observations. 
Therefore, despite the possible contribution from a molecular torus, or other near-neutral plasma 
with low velocity, we still conclude that either the geometry of the X-ray source and/or the 
material which produces the broad, skewed line, differs from source to source. 

The one conclusion which does depend critically on the presence or absence of a narrow line 
is that regarding the rest energy of the disk line. We showed above that a near-neutral origin is 
considerably more likely than a helium-like line if the entire emission line arises from the disk. 
In the presence of a narrow line, we can no longer constrain the energy. We find the fits with a 
6.4 keV disk line are not significantly better than with a helium-like line, if we allow a narrow line 
of typical equivalent width ~ 40 eV at 6.4 keV. This is unsurprising, as the main constraint on 
the rest energy of the disk line comes from the energy of the peak, which in this case is modeled 
with the narrow line. However, given the lack of conclusive evidence for a narrow line component 
in addition to the disk line, a rest energy of 6.4 keV must still be considered more likely and we 
have assumed such an energy in the fits described below. 

4-1-3. Effects of the reflection continuum 

We next explored the effects of Compton reflection on the spectral parameters. As mentioned 
in the introduction, the reflection continuum is commonly observed in the Ginga spectra (NP94) 
and is expected if the material which produces the emission line is optically thick (Lightman & 
White 1988; George & Fabian 1991; Matt et al. 1991). This is is almost certainly the case, given 
the high equivalent widths. 
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Although the reflection component is expected to be relatively weak in the ASCA band, it 
may have some influence on the derived line and continuum parameters (e.g., Mushotzky et al. 
1995; Weaver et al. 1995; Cappi et al. 1996). We therefore added a reflection component to the 
power-law and disk-line, to investigate the possible effects on these parameters. We employed a 
model assuming the abundances and cross-sections of Morrison & McCammon (1983) and elastic 
scattering, following Basko (1978). The fact that the ASCA spectra extend only to ~ 10 keV 
makes the latter approximation reasonable, and computationally expedient. We assumed the 
reflection spectrum from a neutral slab subtending a solid angle of 2n steradian at an X-ray 
source located above the slab. Fits allowing the normalization of the reflection component to be 
free confirmed our suspicions that typically the reflection continuum can be neither detected nor 
constrained for these low-redshift AGN, even when all four instruments are employed. However, 
our hypothesis is that the line arises by fluorescence in an optically-thick accretion disk. Simple 
atomic physics then dictates that the fluorescence must be accompanied by a Compton scattered 
continuum. 

We have therefore repeated our fits with the disk line, including a fixed reflection component, 
assuming an inclination as derived from the disk line model, and an isotropic X-ray source. Note 
that these fits have no additional free parameters. Including such a component can, in principle, 
allow additional constraint to be placed on the disk inclination and may affect the disk-line 
parameters. However, in practice the constraints were no better than without the reflection 
component included. The mean values of q, i and Wk (Table were all consistent in these 
fits compared to those when a pure power-law continuum was used. We did find that differences 
in the line equivalent width tended to be systematically lower, however, with a mean value of 
230 ± 60 eV with reflection included (see Table ||; c.f. 290 ± 50 eV without reflection). The 
reflection component also affects our determination of the continuum (see below). Including the 
reflection component in the fit resulted in a reduction of the total x 2 (i- e - f° r an sources) of 43.9, 
which implies that for the sample as a whole, the reflection model is considerably more likely than 
a simple power-law (Table ^). 

4.2. Kerr model 

The previous section has demonstrated that relativistic effects in an accretion disk orbiting 
a Schwarzschild black hole provide a most agreeable explanation for the line profiles. Those fits 
suggest that the emission is concentrated close to the central black hole. In such an extreme 
environment the details of the black-hole metric may even become important. Laor (1991) has 
presented line profiles for a maximally-rotating black hole in the Kerr metric. It was indeed 
demonstrated that substantial differences in the line profile result, due primarily to the fact that 
the innermost stable orbit is closer to the central point mass (1.23 R g ) and hence the gravitational 
and Doppler effects more intense. 

We have tested such a model on our data including the reflection component, as described 
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above. This time, we have fixed the inner radius at the last stable orbit for a Kerr hole. We fix 
the outer radius at 400 R g , the maximum value allowed in the current implementation of this 
model in XSPEC. The parameters are tabulated in Table [/]. The x 2 f° r these fits (Table ||) are 
remarkably similar to those obtained with the Schwarzschild model given in Table [| The total x 2 
values imply that the Schwarzschild model is ~ 3 times more likely than the Kerr model, but we 
do not consider this to be statistically significant. The similarity is probably a consequence of the 
fact that our fits prefer a relatively low inclination angle. The best-fitting line profiles assuming 
each metric are shown in Fig. ^ for 2 sources with well defined profiles. They are indeed rather 
similar, but the Kerr models have a more pronounced red wing, due to the increased gravitational 
redshift within 6 R g . The mean parameter values for the Kerr fits are entirely consistent with 
those obtained for the Schwarzschild model (Table ^). 

We conclude that a Kerr black hole is not required by our data. Higher-resolution and 
signal-to-noise ratio are required to determine the spin of the black hole definitively. 



4.3. The X-ray continuum 

Whilst our primary goal in considering the 3-10 keV spectra of these Seyfert galaxies was to 
define the iron-line properties of the sample, these data also allow us to define the continuum in 
this band. For a number of the the emission line models mentioned above, we have tabulated the 
mean power-law indices, which are shown in Table ||. As can be seen from this Table, the mean 
power law slope is not strongly dependent on the line model, ranging from < T3_io >= 1.75 for 
a narrow line, to < T3_io >= 1.79 for a disk-line model. However, although the presence of a 
reflection continuum causes minimal change in the derived line parameters there is a substantial 
effect on the continuum. In the fits tabulated in Table |(| the mean spectral slope when reflection 
is included is < T3_io >= 1.91 ± 0.10 (Fig. 1C; Table ||), Ar = 0.12 steeper than with a power-law 
continuum. 

Conceivably there might also be other spectral components which can effect the continuum. 
The most obvious of these are low energy absorption and iron K-edge absorption, both of which 
have been reported for some Ginga spectra (NP94). We have tested for the presence of each of 
these in our sample. Aside from the well-known case of NGC 4151 (for which we have already 
included an absorption column), we find no evidence for soft X-ray absorption. This justifies post 
facto our choice of the 3-10 keV range for determination of the line parameters (note that neutral 
columns <J 10 22 cm~ 2 have little effect above 3 keV). Addition of an iron K-edge does provide a 
significant improvement (at ~ 95 per cent confidence) in 4 observations of 3 sources (NGC 3516 - 
see Fig. |; NGC 3783(1); NGC 3783(2) and NGC 4151(2)). A plausible origin for the edge is in 
the warm absorber mentioned in the introduction. However even when such an edge is detected, 
it has minimal effect on the line and continuum parameters and thus can be disregarded for the 
purposes of our analysis. 
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5. Comparison with previous results 

Although a number of the observations presented here have been published previously, 
detailed discussion of the emission-line properties has been relatively rare. However, where 
comparisons can be made we find consistency with our fits. Specifically, our results are entirely 
consistent with those presented by Ptak et al. for NGC 3227, Mihara et al. for NGC 4051, Weaver 
et al. (1994) and Yaqoob et al. (1995) for NGC 4151, Reynolds et al. (1995b) and Tanaka et al. 
(1995) for MCG-6-30-15, Mushotzky et al. for IC4329A and NGC 5548, Guainazzi et al. for NGC 
7469 and MCG-2-58-22 (Weaver et al. 1995). We find some disagreement in interpretation for 
NGC 3783, where George et al. (1995) claimed no evidence for broadening of the line. We consider 
this to be due to the fact that George et al. fixed the photon index in the fit. Similarly, Leighly et 
al. (1996) claim no evidence for broadening of the line in Mrk 766, but did not allow the energy of 
the line to vary from 6.4 keV. Fig ^ indicates both narrow and broad components to the line in 
Mrk 766. We suggest that Leighly et al. were measuring only the narrow component. 

The results for the ASCA sample are also remarkably consistent with the lower-resolution 
Ginga data. NP94 found a mean energy of 6.37 ± 0.07 keV and equivalent width 140 ± 20 for a 
"narrow" gaussian. As mentioned above, the latter is most directly comparable to that of our 
broad line fits, given the width of the Ginga response, and is entirely consistent with those fits. 
The mean spectral index from Ginga was 1.95 ± 0.05 when the effects of reflection and the warm 
absorber were accounted for (NP94). 



6. Summary of results 

In this paper, we have found: 

• emission at iron Ka is detected in all 23 observations of our 18 sources. 

• the line to be resolved in 80 per cent of the sources, with a mean width corresponding to a 
FWHM of 47,000 km s _1 assuming a gaussian profile. 

• the profiles are often asymmetric, with a peak at ~ 6.4 keV in rest-frame and a broad wing 
extended to the red. The FWZI of the mean line profile corresponds to a velocity of ~ 0.3c. 

• all of the profiles are consistent with that expected from an accretion disk orbiting a central 
black hole. Both Schwarzschild and Kerr metrics model the data equally well. 

• a mean inclination of i ~ 30°. With a random distribution we expect a mean of 60 °. 

• an emissivity function of R~ q between = 6 R g and R a = 1000 R g , with q ~ 2.5. This 
demonstrates that the line emission is concentrated in the inner disk. However, the data are 
not consistent with a universal emissivity law. 
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• assuming the whole of the line arises from the disk, fits with a rest-energy of 6.7 keV are 
strongly disfavored. It is most likely that the line arises by fluorescence in near-neutral 
material. 

• no strong evidence for a narrow component at 6.4 keV in addition to the disk line, but with 
upper limits typically ~ 100 eV. The presence or absence of the narrow component affects 
neither the values of q and i, nor the inferred differences in q, but removes our constraints 
on the rest energy of the disk line. 

• a Compton- "reflection" component, which is expected to accompany the line, is likely to be 
present, but does not affect the inferred line profiles. 

• a mean equivalent width of 230 eV when reflection is included 

• the underlying continuum in our sample to be characterized by a power law with a photon 
index T = 1.9, with a dispersion of 0.15. 

7. Discussion 

Since the discovery of active galactic nuclei, most workers have assumed that they are powered 
by accretion onto a black hole. Such a model apparently presents the only reasonable explanation 
for the observed phenomena (Rees 1984 and references therein). However, as a black hole is not 
directly observable, one is only able to infer its existence through the influence of its physical 
properties: mass and spin. The velocities of gas observed at ~ 0.1 pc from the nucleus of some 
galaxies strongly indicate the presence of huge masses in the central regions (e.g. Ford et al. 
1994; Miyoshi et al. 1995), with the only viable model appearing to be a black hole. However, 
detection of the characteristic relativistic effects of the black hole requires data which probe much 
smaller scales, closer to the event horizon. The velocities here should begin to approach the 
speed of light and, furthermore, the gravitational field of the black hole will cause time-dilation 
effects resulting in gravitational redshift. In order to measure such phenomena, discrete spectral 
features are necessary. The characteristic broadening and distortion of those features due to the 
relativistic motions represent arguably the most compelling test of the black hole paradigm. The 
iron Ka profile observed by Tanaka et al. (1995) in MCG-6-30-15 represents strong evidence for 
such a hypothesis in that object. The similarity of the profile in NGC 4151 (Yaqoob et al. 1995) 
tempts one to extrapolate this model to Seyfert galaxies as a class. But is such a generalization 
appropriate? 

In this paper, we have presented a uniform analysis of a sample of Seyfert 1 galaxies including 
MCG-6-30-15 and NGC 4151. We found results consistent with previous work for MCG-6-30-15 
and NGC 4151, but have further demonstrated that asymmetric line profiles are widespread. Fits 
to a specific physical model for the central regions of the AGN, that of an accretion disk orbiting in 
the Schwarzschild metric, provide significant improvements over a simple gaussian. This supports 



-18- 



the assertion that the line is produced very close to a central black hole. The average emissivity 
function implies that ~50 per cent of the line emission originates within ~20 R g , and ~80 per cent 
within ~100 R g . Special and general relativistic effects thus afford a most satisfying explanation 
for the observed broadening. Furthermore, the lack of any strong blue shifts in our sample argues 
for a planar geometry, such as a disk (e.g., Fabian et al. 1995). However, detailed calculations of 
the line profiles expected in other geometries are beyond the scope of this paper. 

Detailed interpretation of the line profiles is model dependent. However, the distribution 
of parameters for the sample can be compared to those predicted by specific disk models. For 
example, if these Seyfert 1 galaxies were viewed at all inclinations, we would expect a mean of 60° 
for the sample. Therefore, our disk line fits imply that we tend to observe the Seyfert Is at a more 
face-on inclination, on average at 30 degrees. This may be reconciled with Seyfert 1/2 unification 
schemes (e.g. Antonucci & Miller 1985). In these models, all Seyfert galaxies are surrounded by a 
thick, molecular torus which obscures the line of sight in Seyfert 2 galaxies. If the torus axis were 
aligned with that of the disk, we would not then expect to see any edge-on Seyfert 1 galaxies. 
The broad lines would be obscured by the torus for lines of sight with i > i op , the opening angle 
of the torus. Based on our mean disk inclination, we would estimate cosi op ~ (1 — cos < i >/2), 
implying an opening angle of ~ 45°. This is in excellent agreement with estimates based on the 
observations of ionization cones (Pogge 1989 and references therein) in the optical. Although our 
data did not show strong evidence for a narrow emission line - expected from the torus - neither 
did they exclude a narrow component with the predicted strength (Ghisselini et al. 1993; Krolik 
et al. 1993). Conclusive evidence requires data with higher spectral resolution and signal-to- noise 
ratio. However, the presence of a narrow component has little effect on the derived parameters for 
the disk line, particularly q and i. An alternative explanation for the restricted range in i is that 
the disk itself - if geometrically thick - obscures the edge-on lines-of-sight. 

Assuming that the disk line accounts for all of the emission at iron Ka, the requirement 
for a rest-energy of 6.4 keV implies that the line arises by fluorescence excited by the power law 
continuum. This raises the question as to how the disk can remain relatively neutral when exposed 
to such intense illumination. The immediate conclusion is that the density of the fluorescing 
material is extremely high. Photoionized accretion disks have been discussed by Ross & Fabian 
(1993) and Matt, Fabian & Ross (1993). For a standard a-disc, we expect a density, n, at 20i? g , 
of (Matt et al. 1993): 

2 

n = 1.2 x 10 15 ■ - %06 cm- 3 

M^Mgaom 

where M= O.lMo.i is the accretion rate in units of the Eddington rate, M = 10 8 M§ is the 
black hole mass in solar masses, rj = 0.06r/o.06 is the energy conversion efficiency and a = O.Olao.oi 
is the dimensionless viscosity parameter. The ionization parameter is approximately: 
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^ ~ nR 2 ~ M 8 r]l 06 

where Lx = 10 44 L44 is the ionizing luminosity. Below £ ~ 300 a rest energy close to 6.4 
keV is predicted for the line (Matt et al. 1993). Thus even in the central regions, we typically 
expect the mean ionization state for the disk to be low. However, presumably the surface layers, 
where the bulk of the fluorescence occurs, could be more highly ionized than is indicated by the 
above relations. Assuming fixed energy fractions between the X-ray and intrinsic disk emission, 
we infer that £ ocM 3 and thus that the accretion rate is the dominant parameter in determining 
the ionization state of the disk. This implies that the M values in these Seyfert galaxies are 
sub-critical. Interestingly, for some quasars, in which the accretion rates may be higher, there is 
evidence that the peak of the emission line is higher, indicative of more ionized iron (Nandra et 
al. 1996a). 

Our fits to the Schwarzschild model have also shown that the measured line profiles are 
not compatible with a universal emissivity function, even when the inclination and a possible 
contribution of a narrow line from another region are allowed to vary. This might be a consequence 
of slight differences in geometry from source to source. For example, if the true inner and outer 
radii of the disk are different to those assumed, this can be compensated for in the fit by changing 
the apparent emissivity function. Alternatively, as we are dealing with a fluorescence line, the line 
emissivity depends upon the X-ray illumination, which could change from source to source. For 
example, with a point-like X-ray source located above the disk, the line profile depends upon the 
height of the source (Matt et al. 1992). Our data are only well- fit by the disk line model if we allow 
the emission to be distributed over a range of radii. Emission at a single radius cannot account 
for the observed profiles in the majority of cases. This argues against models where the height of 
the X-ray source is very small, where only the inner disk is illuminated. Furthermore, models in 
which the disk presents significant solid angle to the X-ray source only at a small number of radii 
are difficult to reconcile with our data. One might envisage such a scenario if the X-ray source 
were at the inner edge of the disk itself or if the very innermost regions of the disk were puffed up 
by radiation pressure. 

A simplistic "one zone" model for the X-ray source may not be appropriate. A currently- 
popular model is of multiple regions or "hot spots" covering the surface of the disk (e.g., Stern 
et al. 1995). The mean emissivity profile is harder to calculate in such a case, as it depends 
on the distribution of hot spots within the X-ray producing region and the variation of their 
density with radius. By observation, we conclude that in such a model, this latter quantity should 
vary approximately as i?~ 2 ' 5 . The challenge is now for theoretical models to reproduce such a 
relationship and, one would hope, relate it to the energy-generation mechanism for the X-rays. 

An area where our observational data come into mild conflict with the conventional wisdom 
regarding disk reflection is that of the line equivalent width. The detailed calculations of George 
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& Fabian (1991) and Matt et al. (1991) predicted a value from the disk of Wk = 140 eV for a 
continuum with T = 1.9 and an inclination of 30°. This is a factor 1.2 — 2 less than our observed 
mean of 230 =L 60 eV. An obvious way to resolve this discrepancy is to revise our assumptions 
about the metal abundances. The required increase in equivalent width can be achieved by an 
enhancement in the iron abundance of a factor 1.3-3.5 over the value of Morisson & McCammon 
(1983) which are usually employed (see Fig. 16 of George & Fabian 1991). Indeed, using the 
updated solar abundances of Anders & Grevesse (1989), Reynolds, Fabian & Inoue (1995a) 
find that equivalent widths of ~ 250 eV could be achieved without any enhancement in iron. 
The possible contribution of ~ 30 eV of a narrow line from another region would also ease this 
discrepancy somewhat. 

The mean iron abundance derived here argues that the material being accreted is relatively 
metal-rich, consistent with recent broad-line models for quasars (Hamann & Ferland 1993). Since 
the central regions of most spiral galaxies have similar metallicities, but the abundances drop 
with increasing radius (McCall, Rybski & Shields 1985; Shields, Skillman & Kennicutt 1991). 
This tends to imply that the accretion fuel is produced by stellar mass-loss relatively close to 
the nucleus and is not brought in from further out in the galaxy. The required accretion rate 
is M= 3 x 10~ 2 L447/o.o6 M yr -1 (here L44 represents the accretion luminosity, rather than the 
ionizing luminosity). The mass- loss rate is expected to be M~ 10 _2 Mio MQyr -1 , where M\q is 
the mass of stars in units of 10 10 M , which seems plausible. 

Our fits have also allowed us to define the continuum shape in the 3-10 keV band. We find 
a mean slope of ~ 1.8 with a power-law plus line model. This lies between the values determined 
by HEAO-1 and EXOSAT (using a single power-law; Mushotzky 1984; Turner & Pounds 1989) 
and that determined from Ginga, where a more complex model was employed, including a warm 
absorber and reflection. When we include the reflection component in our ASCA fits, we find 
r = 1.9, in excellent agreement with the Ginga data, which were arguably more appropriate for 
continuum determination. These results indicate that the reflection component has a substantial 
effect on the determination of the continuum using ASCA data in the 3-10 keV band, but that 
the warm absorber does not. We will return to this latter point in paper III. Our data have 
implications for the physical models of the X-ray continuum in AGN. For example, the disk/corona 
model of Haardt & Maraschi (1993) has been successful at explaining the mean intrinsic slope for 
Seyfert galaxies of V = 1.9. However, Stern et al. (1995) have pointed out that this model has 
some difficulties in explaining a range of observed slopes, and suggest a solution in which a number 
of distinct active regions illuminate the disk. Our data confirm the range of slopes observed by 
Ginga, and thus lend support to this hypothesis. 
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Table 1. The ASCA Seyfert 1 sample. 



Name 


RA 


DEC 


Redshift 


A H (Gal) a 


N obs/ N tot 


logZ& 




(1950) 


(1950) 


(V93) 


10 20 cm- 2 




erg s -1 


Mrk 335 


00 03 45.1 


19 55 29 


0.025 


4.0 


1/1 


43.42 


Fairall 9 


01 21 51.0 


-59 03 54 


0.046 


3.0 


1/1 


44.26 


3C 120 


04 30 31.5 


05 14 59 


0.033 


12.3 


1/1 


44.34 


NGC 3227 


10 20 46.9 


20 07 06 


0.003 


2.2 


1/1 


42.01 


NGC 3516 


11 03 22.8 


72 50 18 


0.009 


2.9 


1/1 


43.43 


NGC 3783 


11 36 33.0 


-37 27 42 


0.009 


8.9 


2/2 


43.25 


NGC 4051 


12 00 36.3 


44 48 35 


0.002 


1.3 


1/1 


41.56 


NGC 4151 


12 08 00.9 


39 40 59 


0.003 


2.1 


3/6 


42.97 


Mrk 766 


12 15 55.7 


30 05 25 


0.012 


1.6 


1/1 


43.08 


NGC 4593 


12 37 04.6 


-05 04 11 


0.009 


2.0 


1/1 


43.06 


MCG-6-30-15 


13 33 01.4 


-34 02 27 


0.008 


4.1 


2/2 


43.07 


IC 4329A 


13 46 28.2 


-30 03 42 


0.016 


10.4 


1/1 


43.94 


NGC 5548 


14 15 43.6 


25 21 58 


0.017 


1.7 


1/1 


43.76 


Mrk 841 


15 01 36.2 


10 37 53 


0.036 


2.2 


2/2 


43.82 


NGC 6814 


19 39 55.7 


-10 26 35 


0.006 


9.8 


1/3 


41.28 


Mrk 509 


20 41 26.1 


-10 54 16 


0.035 


4.2 


1/1 


44.38 


NGC 7469 


23 00 44.4 


08 36 17 


0.017 


4.8 


1/3 


43.60 


MCG-2-58-22 


23 02 07.2 


-08 57 20 


0.047 


3.4 


1/1 


44.14 



a Galactic Hi column density from 21cm measurements 

b Number of observations for which adequate data were obtained/Total number of observations 
c Mean 2-10 keV luminosity 

References. — V93: Veron-Cetty & Veron (1993); Ah values are from Elvis, Lockman & Wilkes 
(1989); Stark et al. (1992) or the HEASARC online service. 



-23- 



Table 2. Narrow line fits 



Name 


A a 


1 3-10 


(keV) 


rci 

J Ka 


W Ka 

(eV) 


X NL/d.O.f. 


El/ 

^2 


Mrk 335 


0.29 


1.87 ± 0.10 


6 37+2-Jf 

uo ' —0.16 


1 2 + 0"I 

x ' -0.8 


110+go 

u -70 


227.9/252 


5.0 


Fairall Q 


60 


1 oo i n nc 

1 .UU _L_ \J . \ ' \J 


6 3S+ - 15 

u -°°-0.08 


9 Q+ - 9 
z - y -1.0 


120+fS 

J-^ u -40 


640 6/610 


15 


3C 120 


1.31 


1 79 ± 03 


uo ' -0.05 


3 h +0 A 

°- tJ -1.7 


,u -30 


881 9/950 


15.9 


NGC 3227 


0.48 


1.52 ± 0.04 


6 42 +0 ,-°j 


o c+0.9 
o,tJ -0.9 


120+12 


885 3/920 


30.2 


NGC 3516 


2.03 


1.74 ±0.03 


f- O7+0.03 
u -°'-0.03 


7.9+fj 


90l?8 


1055.3/1004 


42.8 


NGC 3783(1) 


1.03 


1.57 ±0.05 


c qo+0.04 
b - ri8 -0.04 




150+1° 


800.8/788 


34.2 


NGC 3783(2) 


1.11 


1.51 ±0.05 


6 16+ - 14 


5.71U 


80118 


715.1/743 


27.4 


NGC 4051 


0.73 


1.92 ±0.06 


6 44+ - 07 
°-^ 4 -0.07 


2.8l?; 2 7 


1401^8 


672.4/657 


13.3 


NGC 4151(2) 


3.75 


1.42 ±0.12 


6 40+ - 02 
D -^ u -0.03 


29+^ 


1 1 +2 ° 


760.4/693 


63.2 


NGC 4151(4) 


6.31 


1.59 ±0.07 


6 S6+ - 05 


32±§ 


90+ 30 

yu -20 


1409.1/1401 


47.2 


NGC 4151(5) 


7.73 


1.65 ±0.05 


c o 7 +0.02 
D -°' -0.02 


3111 


90+ 20 

yu -20 


1348.6/1228 


45.0 


Mrk 766 


0.75 


2.00 ±0.05 


c O7+0.03 
D " 3 ' -0.03 


z -°-0.7 


140118 


584.3/682 


21.2 


NGC 4593 


0.92 


1.78 ±0.05 


6-3518:81 






827.4/868 


12.7 


MCG-6-30-15(l) 


1.64 


1.94 ±0.05 


6.26i8;? 7 


3.2il:| 


70±|o 


943.4/845 


9.0 


MCG-6-30-15(2) 


1.11 


1.81 ±0.05 


6.38t8:? 7 


a O+0.9 
d.»_ 20 


iool 2 8 


829.7/802 


9.7 


IC 4329A 


1.95 


1.70 ±0.03 


I? on+0.05 
°- oz -0.05 


7 3 +L4 

' -°-1.6 


80l 2 8 


1297.7/1249 


38.5 


NGC 5548 


1.24 


1.75 ±0.04 


fi on+0.04 
u - oy -0.04 


4 6 +L3 


qn+30 
yu -30 


973.6/1000 


21.6 


Mrk 841(1) 


0.38 


1.84 ±0.09 


6 37+ - 12 

u -°'-0.02 


2 4 +o.9 

z -^-0.8 


170I88 


378.1/368 


7.0 


Mrk 841(2) 


0.21 


1.57 ±0.10 




1-41^ 


12018*8 


181.9/203 


4.7 


NGC 6814(1) 


0.02 


1.54 ±0.26 


5.9+™ 


0.418J 


230ll?8 


68.1/76 


4.1 
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Table 2 — Continued 



Name 


A a 


1 3-10 


(keV) 


jd 


(eV) 


X 2 NL/d.o.f. 


H 


Mrk 509 
NGC 7469(2) 
MCG-2-58-22 


1.21 
0.76 
0.29 


1.76 ±0.05 
1.78 ±0.07 
1.58 ±0.08 


6 S7+ - 07 

u -°'-0.07 
a oq+0.06 

o-^-o.o? 


o r:+0.7 
o.O_ 2 4 

3.811:1 

2.61SJ 


13011° 
150±S8° 


819.0/896 
501.1/451 
408.4/388 


7.4 
19.7 
7.8 



a Power law flux at 1 keV; 10 3 ph cm 2 s 1 keV 1 
b Power-law index 

c Energy of the emission line in the rest-frame of the source 
d Flux of the line; 10 -5 ph cm~ 2 s _1 
e Line equivalent width 

f F-statistic for the addition of two free parameters 

Note. — Fits were undertaken in the 3-10 keV range. Absorption is by Galactic N-r only, except 
in the case of NGC 4151, where the column was left free. Error bars on the fit parameters are 68 
per cent confidence limits for 3 interesting parameters (Ax 2 =3.5). 
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Table 3. Broad line fits 



Name 


1^3-10 










ivirK ooD 








< 1.1 


n q 
u.y 


Fairall 9 


i n-i +0 08 




or,n + 130 

320l^ u 


„ or +0 18 


21.7 




i on+0.08 
i -° y -0.06 


« on+0.20 

6.39I 020 


oon+200 
•""-ISO 


0.74_ 027 


on n 


i\Tri /~i ooo'v 


1 co+0.08 
i - oz -0.03 


6.30_ 012 


i on+60 
180_ 50 


n oc+0.12 
u - zo -0.10 


7.0 


MP p o cr i 

INLtU oolo 


-i oq+0.05 


fi i 9+0.13 


d/in+130 


n R/1+0-22 
°- 64 -0.19 


/in i 
4U.1 




-I 7n +0.11 

L70 -o.o8 


fi 1 1+0-17 
6 - n -0.18 


510 -190 


n 71 +0.35 
0-71-0.28 


i n q 


1VPP ^7^(9\ 




r nq+0.26 
d - yo -0.39 


q fin +300 


n ?n +0.48 
U -' U -0.36 


9/1 ^ 


JNCjtU 4U51 








< 1.4 


1 o 

1.8 


4151(2) 


i ^+0.14 
1- 44 _0.13 


a on+0.05 


i on+60 

130_ 40 


n 1 n+0.08 

O-10-o.io 


4.3 


INGC 4151(4) 


1 crv+O-lO 

l- 57 lo.n 


a 1 ,1+0.13 
6 - 14 -0.18 


Q/in+160 


n cc+0.25 
0- 56 -0.19 


cor 1 

52.6 


IN^jrL/ 4101^0 J 


1.58_ n 


cr qo+0.12 
°- yz -0.17 


qqn+140 


n «/| +0.22 

U.b4_ 15 


DU.l 


Mrk 766 


— U. 11 


6.03+8|i 


550 oQn 


0.72+8I? 

— U.zo 


16.0 


NGC 4593 








< 1.1 


3.4 


MCG-6-30-15(l) 


2-o4±g:S? 


6.06l8f 


320l&° 


o.62!8:^ 


27.3 


MCG-6-30-15(2) 


1 - yo -0.08 


5-80±8;» 


4601*8 


o.88i8;l 


32.1 


IC 4329A 


i 71 +0.03 
11 ' J'-O.OS 




nn +40 
iiu -30 


u - 14l -0.08 


10.7 


NGC 5548 


1 7 q+0.05 
L - ' y -0.05 


6 41 +0 - 16 
D - 4i -0.14 


17n +60 
LIU -40 


u -°^-0.12 


5.8 


Mrk 841(1) 


z - uu -0.15 


6-is±8:li 




n fi o+0.35 
U - Od -0.32 


6.9 


Mrk 841(2) 








< OO 


0.9 


NGC 6814(1) 


1 q+2.0 


6.318:1 




o.6i8i 


8.9 


Mrk 509 


-i on+0.09 
L -° z -0.07 


6.351?$ 


2101118 


o.4T+8:i8 


14.2 


NGC 7469(2) 




6.39l°;It 


2801$ 


S7 +0 - 31 
u -°' -0.37 


3.7 


MCG-2-58-22 








< 0.85 


1.7 



a Width of the gaussian (keV) 

b F statistic for the addition of a as a free parameter 

Note. - Fits with a power-law plus broad-line to the ASCA data in the 3-10 keV range. 
Absorption is by Galactic Nu only, except in the case of NGC 4151, where the column was left 
free. Error bars on the fit parameters are 68 per cent confidence limits for 4 interesting parameters 
(A x 2 =4.7). 
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Table 4. 


Disk line fits 


; Schwarzschild geometry; 


free emissivity 




1 > cLIIlfc! 


1 3-10 






Q 




Mrk 335 




oon+280 
ZJU -170 


oo+67 
ZJ -23 


2- z+ S 

— OO 


1.3 


Fairall 9 


— U.Uo 


ocn+180 
oou -130 


451M 


1.9l°; 5 7 


-0.8 


3C 120 


1.87t°-«j 

— U.UO 


290i 7 J° 


60l?° 


9 9+0^3 
5 


-3.0 


NGC 3227 


—0 05 


210+ 80 
ziu -80 


2011° 


2-Olgl 

— u.o 


4.4 


NGC 3516 


i.80+g : gt 




27l£ 




21.5 


NGC 3783(1) 


l-64±8-g§ 

— u.uu 


o 7^+330 
o/u_ 110 


2ll 7 ° 


2.5+~ 4 


10.2 


NGC 3783(2) 


1 co+0.07 

— U.U f 


430+ 140 


09+3 

az -13 


5.2+?° 8 

— 1.0 


8.0 


NGC 4051 


2.00+Hs 

— U.Uo 


46011° 


33l? 3 


4.4+^ 


13.0 


NGC 4151(2) 


— U.lO 


190l|° 


1811 


2.0+^ 


3.9 


NGC 4151(4) 




oou -120 


3311 


5.oll:l 


19.8 


NGC 4151(5) 


— u.uy 


3201™ 


17l| 


Z.o n q 
— U.o 


28.5 


Mrk 766 


o 10+O.II 
— U.Uo 


6101158 


3511 


3.7+1- 9 

— l.U 


7.9 


NGC 4593 


-. OO+0.08 
J--00 q Qy 


24n +i40 

Z4lu -140 


45+ 45 
40 -45 


2.21^ 


-0.3 


MCG-6-30-15(l) 


2.041°^ 


430+ 160 

4OU -160 


33^ 


1+1.1 

—0 6 


7.8 


MCG-6-30-15(2) 


1 qi +0.08 
1 - yi -0.06 




34l5 


4.2+f° 


3.0 


IC 4329A 


i 7 r,+0.03 
L - ' z -0.03 


1501*° 


99+15 
zz -22 


2-ilS:S 


2.6 


NGC 5548 


1 7Q+0-05 
1 -' y -0.05 


200+ 9 7 ° 


47l£ 


1.9I8:? 


-0.6 


Mrk 841(1) 


i qc;+0.l4 


51011° 


271? 


2.6l|° 


6.4 


Mrk 841(2) 


1 62+ - 16 


500115° 


38l? 2 


1Q +oo 


3.6 


NGC 6814(1) 












Mrk 509 


1 81+0- 11 


lTOlJg 


41 -32 


2-i±S:S 


-1.4 


NGC 7469(2) 


i-8o±8:S§ 


aiolgg 


iq+71 
iy -19 


1-81^° 


-1.8 


MCG-2-58-22 


1.601H? 


160ir 


71 + 19 
IL -71 


-iol^ 


0.3 



inclination of the disk to the line-of-sight 
b The emissivity varies as R~ q 

c Difference in x 2 compared to the broad-line fit (x 2 bl~ X 2 s) 

Note. — Absorption is by Galactic Nn only, except in the case of NGC 4151, where the column 
was left free. Error bars on the fit parameters are 68 per cent confidence limits for 4 interesting 
parameters (Ax 2 =4.7). 
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Table 5. Disk line fits; Schwarzschild geometry; q=2.5 



Name 


1^3-10 




i 


pa 


X 2 S25/d.o.f 


IVlrK OoO 


i qo+0.14 


9 o n +160 


9zL +13 
24 -24 


n 9 


997 Q /9^9 

LL i .y / zoz 


ran an y 


1 qi+0-08 


oqn+130 


09+12 


0.0 


fi9^ 3 /fii n 
ozo.o / oiu 


3P 1 9fl 

Ov 1ZU 


1 fiO+0-05 
i - 8y -0.04 


otrn+90 
JOU -100 


59+ 10 
oy -13 


0.0 


8^9 7/0.^0 


wpp 9.007 
IN LiL/ o/Z f 


L5b -o.o5 


260_ 60 


90+8 


0.0 


Q70 /oon 
(y.y/yzU 


INvjrLy OOIO 


i - 8U -0.04 




9fi+ 4 


1 A 
1.4 


QQ/1 n /1 nn/i 
yy4.u/ iuu4 


WPP ^S^f 1 ^ 


l.O4_ 005 


O 7n + 70 


91+5 


n n 
u.u 


771 A /788 
( 1 1.4/ ( 00 


npp "^ys^ 


i cc+0.06 
i - oo -0.05 


otrn+60 


q+H 

y -9 


1 n fi 


RQ4 3 /743 


wpp /in^i 

IN^j^ 4U01 


i n 7 +0.06 
J- 9 ' -0.06 




97+7 


1 1 
0. 1 


oou.y / 00 1 


WPP 41 ^1 


i - 4i -0.12 


ocn+60 
2DU_ 50 


9H+5 

—5 


0.0 


7^8 9 /fiQ4 

i oo.z/ oy4 


WPP 41 ^1 f4"\ 

IN^JIV^ 410114 1 


i co+0.08 
i - oz -0.07 


9on+ 50 

ZyU -40 


91+5 


0. 1 


1 343 s /1 4nn 

1040.0/ 14UU 


NPP 41 M fti 


i 6n+ 07 

i - DU -0.07 


250+ 40 


15+ 4 


7 4 


1 97^ 1 /I 997 


Mrk 766 


2-0918:8? 


440l{^ 


35i 


9.2 


570.7/682 


NGC 4593 


1-8418:88 


280^8 


4511? 


1.1 


825.5/868 


MCG-6-30-15(l) 


9 n9 +0.05 
z - uZ -0.05 


o 2f) +70 
JZu -100 


34±| 


4.5 


910.8/845 


MCG-6-30-15(2) 


-i 07+0. 06 
J-- 8 ' -0.05 


2601*8° 


331^ 


14.4 


809.2/802 


IC 4329A 


1 79 +0.03 
x - ' z -0.02 


180l|8 


26l 8 7 


4.1 


1289.1/1249 


NGC 5548 


1 7O+0.05 
L - ' y -0.04 


220188 


3911? 


4.5 


973.0/1000 


Mrk 841(1) 


1 qr;+0.l0 


5001^8 


2718 


0.1 


364.8/368 


Mrk 841(2) 


1 61 + - 12 

^^-O.ll 


2501188 


30l? 5 


4.4 


181.4/203 


NGC 6814(1) 












Mrk 509 


1 01 +0.06 
i - 8i -0.05 


210U88 


401^ 


2.0 


809.4/894 


NGC 7469(2) 


i - 8b -0.09 


370+iso 


451^ 


1.7 


500.7/451 


MCG-2-58-22 


1 62+ - 11 
^•"^-o.os 


0-1 n+ 170 


4111° 


0.5 


406.8/388 



a F statistic for the addition of q as a free parameter 

Note. — Absorption is by Galactic iVn only, except in the case of NGC 4151, where the column 
was left free. Error bars on the fit parameters are 68 per cent confidence limits for 3 interesting 
parameters (Ax 2 =3.5). 
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Table 6. Disk line fits; Schwarzschild geometry; with reflection 



AT 


1 3-10 




i 


9 


A 2a 


Mux 


Mrk 335 


2.05+2^ 

—0.23 


230+? 40 , 

— 190 


99+68 
—22 


2.3+°° 


1.1 


0.95 ± 0.12 


Fairall 9 




320+?f° 

"->^ u — no 


46 +4 i 


1 9+°'? 


0.6 


2 05 ± 27 


3C 120 


J - -CM —0.14 


260+^° 


60+?° 

""—13 


9 O+0.4 
°— 0.6 


0.0 


4.69 ± 0.61 


NGC 3227 


1 67+2-21 

1,u ' -0.05 


190+^2 


20+?, 

14 


2 0+?'£ 

^• u -1.0 


3.9 


2 64 ± 33 


NGC 3516 


1 OQ+0.04 
JU -0 .04 


320+^° 

" — 90 


27+v 

1 -7 


2.8+^°, 

u —0.4 


-2.5 


6.82 ± 1.03 


NGC 3783(1) 


i 77 +0.08 
' ' —0.06 


590+}^ 

"^"—140 


35_-^g 


8.5+?° R 

u • — 4.6 


2.9 


5.12 ± 0.67 


NGC 3783(2) 


1.70+H? 

—0.0 1 


400+112 


32+?. 

— 15 


6.4+2° fi 


2.0 


6.12 ±0.80 


NGC 4051 


9 i 1 +0.08 
^- 11 — 0.08 


420+i?2 


34+?, 


5 1 + S° C 


0.6 


2 13 ± 28 


NGC 4151(2) 


1.60+2"^ 

""—0.13 


150+?° 

""—70 


17+1? 

— 17 


1.8+°f 


5.6 


24.3 ± 3.1 


NGC 4151(4) 


1 54+°-J? 

- L -°'±-0.11 


490+i2° 
^ u -90 


33+?o 


5 4+?-° 

1.5 


2.3 


27 8 ± 3 6 


NGC 4151(5) 


i 7 9+0.09 
1- ' ^—0.08 


290+12 

60 


17+!? 

x ' -5 


2 9+2'? 

0.3 


1.1 


32 6 ± 4 2 


Mrk 766 


9 9C-+0.10 
z, - z,u — 0.09 


580+^22 


34+3 


3 9+?'2 
°- J -1.0 


1.4 


1 95 ± 25 

J- ■ t_/ .J — i — V 7 • — 7 


NGC 4593 


1.91+2-2? 

-0 .05 


90+12 

—50 


—0 


1.5+^ 3 

A '" — oo 


5.1 


3.38 ± 0.44 


MCG-6-30-15(l) 


9 1fi +0.06 
z - 1D -0.07 


4ooiyg 


33+J 


o 4 +1.3 
J - 4 -0.8 


3.2 


4.69 ±0.61 


MCG-6-30-15(2) 


2 ns +0 - 07 

z - uo -0.07 


410ljg 


34+j 


44 +3.o 


-0.4 


3.83 ±0.50 


IC 4329A 


i o 4 +0.03 
i - 84 -0.03 


120+™ 


17 +14 
1<_ 17 


9 1+2^8 
z - 1 -l.l 


3.3 


8.12 ± 1.06 


NGC 5548 


1 89+ - 04 
i - oy -0.05 


8011° 


oi 76 


-9.0™ 


6.2 


4.71 ±0.61 


Mrk 841(1) 


2 n7 +o.i3 

z - u ' -0.13 


470ii° 


27+? 


2.61- 


1.0 


1.29 ±0.17 


Mrk 841(2) 


1 74+0.I6 

x - /4 -0.13 


440+1° 




10+°° 

iu — oo 


-1.7 


1.07 ±0.14 


NGC 6814(1) 












0.13 ±0.02 


Mrk 509 


i 09+O.O6 


110+ 7 °° 


97+63 
z ' -27 


9 9+1.2 
z - z -oo 


1.6 


4.58 ±0.59 


NGC 7469(2) 


1 Qi +0-09 
± - yi -0.08 


120lgg 


n +89 
u -0 


-9 0+ 10 ' 7 


4.1 


2.78 ± 0.36 


MCG-2-58-22 


i 7o+°- 16 

± -' u -0.18 


1401S8 


4D -46 


-4 8+°° 

^• u — oo 


2.0 


1.48 ±0.19 



a Reduction in x 2 relative to the model with a Schwarzschild disk line and power law, i.e. X 2 S-X 2 SR 

b Flux in the 2-10 keV band in units of 10 -11 erg cm -2 s _1 , based on the model described in this 
table, but corrected for absorption. The error bar corresponds to the absolute calibration of the 
SIS detectors of ~ 13 per cent 

Note. — Absorption is by Galactic Nr only, except in the case of NGC 4151, where the column 
was left free. Error bars on the fit parameters are 68 per cent confidence limits for 4 interesting 
parameters (Ax 2 =4. 7). 
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Table 7. Disk line fits; Kerr geometry; with reflection 



Name 


T3-18 






q 


A x 2a 


Mrk 335 


2 04 +0 /^ 


^ Ou -220 


22+S 


9 9+I.2 
z - z -oo 


-0.6 


FViirnll Q 

_L Ojll Ojll u 


1 OQ+0.14 

J^-^-o.io 


o fin +200 
^ DU -150 


89 +1 
oy_ 49 


I q+0.9 

1.0_ 7 2 


n 1 


3P 1 9D 


z - ui -0.04 


fi7n +90 
°' u -120 


88+ 2 


9 7+0.4 
z -'-0.2 


5 3 
0.0 


wpp 3997 


-1 fi7 +0.05 


io n +110 

i»u_ 70 


91+7 
zi -21 


1 8+ - 6 


-l.o 


WPP 




oon+90 
<JO"_90 


96+ 3 


2 6 +0 - 3 

z - D -0.2 


4 Q 


wpp 3783(1 


77 +0.06 
1 - ' '-0.06 


550+ 190 
oou -190 


26 +5 


9 7+0.3 
z -'-0.4 
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a Reduction in x 2 relative to the model with a Schwarzschild disk line and reflection continuum, 
i-e. X 2 sr-X 2 kr 

Note. — Absorption is by Galactic Nr only, except in the case of NGC 4151, where the column 
was left free. Error bars on the fit parameters are 68 per cent confidence limits for 4 interesting 
parameters (Ax 2 =4.7). 
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Table 8. Comparison of \ 2 values 



Name 


X BL 


X S 


2 

X S67 


2 

X ss 


2 

X SR 


2 

X KR 


Q.O.I 


Mrk 335 


227.0 


225.7 


226.3 


225.3 


224.6 


225.2 


251 


Fairall 9 


618.6 


619.4 


625.5 


626.9 


618.8 


618.9 


609 


3C 120 


846.3 


849.3 


868.9 


862.7 


849.3 


844.0 


949 


NGC 3227 


878.0 


873.6 


890.3 


875.9 


869.7 


871.5 


919 


NGC 3516 


1014.1 


992.6 


994.5 


1011.3 


995.1 


990.2 


1003 


NGC 3783(1) 


781.6 


771.4 


775.8 


771.4 


768.5 


770.1 


787 


NGC 3783(2) 


692.4 


684.4 


683.1 


693.2 


682.4 


683.2 


742 


NGC 4051 


670.8 


657.8 


656.4 


645.0 


657.2 


657.0 


656 


NGC 4151(2) 


755.9 


752.1 


763.8 


757.1 


746.5 


749.7 


693 


NGC 4151(4) 


1357.7 


1337.9 


1331.6 


1347.5 


1335.6 


1336.2 


1399 


NGC 4151(5) 


1295.9 


1267.4 


1280.9 


1282.1 


1266.3 


1266.3 


1226 


Mrk 766 


570.9 


563.0 


562.2 


577.2 


563.0 


563.5 


681 


NGC 4593 


824.2 


824.5 


831.7 


824.5 


819.4 


819.9 


867 


MCG-6-30-15(l) 


913.8 


906.0 


905.5 


898.6 


902.8 


903.3 


844 


MCG-6-30-15(2) 


797.7 


794.7 


792.1 


797.9 


795.1 


792.2 


801 


IC 4329A 


1287.5 


1284.9 


1296.2 


1306.9 


1281.6 


1281.1 


1248 


NGC 5548 


968.0 


968.6 


974.2 


974.1 


962.4 


964.1 


999 


Mrk 841(1) 


371.1 


364.7 


369.9 


359.4 


362.3 


363.4 


367 


Mrk 841(2) 


181.1 


177.5 


180.9 


181.8 


179.2 


182.1 


202 


NGC 6814(1) 
















Mrk 509 


806.2 


807.6 


816.7 


809.7 


806.0 


805.2 


893 


NGC 7469(2) 


497.0 


498.8 


503.3 


498.4 


494.7 


494.9 


450 


MCG-2-58-22 


406.6 


406.3 


407.4 


407.8 


404.3 


404.3 


387 



Total Values 



Value (-16000) 762.4 628.2 737.2 734.7 584.4 586.3 973 

Probability 0.87 0.97 0.90 0.90 0.98 0.98 
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Table 8 — Continued 



Name x 2 bl X 2 s X 2 S67 X 2 ss X 2 sr X 2 kr d.o.f 



Relative likelihoods" 

> 10 3 > 10 3 

> 10 3 > 10 3 

> 10 3 > 10 3 

> 10 3 > 10 3 
0.39 

2.6 



a Likelihood of the model in the column heading compared to the row, based on the total x 2 
values. 

Note. - The subscripts used are: BL - broad gaussian (Table |||); S - Schwarzschild model 
(Table |j); S67 - Schwarzschild model with a rest energy of 6.7 keV; SS - Schwarzschild model 
at a single radius; SR - Schwarzschild with reflection (Table |^); KR - Kerr model with reflection 
(Table g) 



X BL 






x s 


< 10" 


-3 


^2 

X S67 


< 10" 


-3 


X 2 ss 


< 10" 


-3 


X 2 SR 


< 10" 


-3 


X KR 


< io- 


-3 



> 10 3 

> 10 3 

> 10 3 

< 10- 3 

< 10~ 3 



> 10 3 

< 10~ 3 

3.5 

< 10- 3 

< 10~ 3 



> 10 3 

< 10" 3 
0.29 

< 10~ 3 

< 10" 3 
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Table 9. Mean parameters for the Seyfert 1 galaxies 



Parameter Unit ji a ji b < fi > c of 



Narrow Line Fits, 3-10 keV 



T3-10 




1.74 ±0.14 


1.74 ± 0.01 


1.75 ±0.05 


0.12 ±0.04 




keV 


6.34 ±0.12 


6.36 ±0.01 


6.36 ±0.02 


< 0.06 


w Ka 


eV 


115 ±39 


98 ±8 


98 ± 12 


< 24 


Broad Line Fits, 3-10 keV 


F3-10 




1.80 ±0.17 


1.76 ± 0.02 


1.78 ± 0.06 


0.13 ±0.05 


ERa 


keV 


6.29 ±0.14 


6.34 ± 0.02 


6.34 ±0.04 


< 0.09 


W Ka 


eV 


290 ± 240 


150 ± 20 


160 ± 30 


< 90 


<?Ka 


keV 


0.40 ± 0.28 


0.39 ± 0.06 


0.43 ±0.12 


< 0.29 


Disk Line, free q 


F3-10 




1.80 ±0.17 


1.76 ± 0.02 


1.79 ± 0.07 


0.15 ±0.05 


w Ka 


eV 


310 ± 140 


280 ± 30 


290 ± 50 


< 140 


i 


o 


35 ± 15 


29 ± 2 


29 ±3 


< 9 


q 




1.9 ±3.2 


2.5 ±0.2 


2.5 ±0.4 


< 0.8 


Disk Line, free q, reflection 


F3-10 




1.92 ±0.17 


1.88 ± 0.02 


1.91 ±0.07 


0.15 ±0.05 


w Ka 


eV 


270 ± 160 


180 ± 20 


230 ± 60 


110 ±50 


i 


o 


26 ± 17 


29 ± 2 


29 ±5 


< 11 


q 




1.1 ±4.5 


2.8 ±0.3 


2.8 ±0.4 


< 0.9 
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Table 9 — Continued 



Parameter Unit fx a Li b < ll > c of 



Kerr model, free q, reflection 



T3-10 




1.90 ±0.17 


1.88 ± 0.02 


1.89 ±0.07 


0.15 ±0.05 


w Ka 


eV 


320 ± 230 


210 ± 30 


270 ± 60 


130 ± 60 


i 


o 


34 ±28 


28 ±3 


28 ±4 


< 8 


Q 




1.3 ±3.1 


2.8 ±0.2 


2.8 ±0.2 


< 0.4 



a Unweighted mean and standard deviation 
b Weighted mean 

Expectation allowing for measurement errors (see text) 
d Intrinsic dispersion of distribution (see text) 

Note. — Error bars are 68 per cent confidence limits. Upper limits are 90 per cent confidence. 
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Fig. 1. — Histograms of the line energy (left panels) and the equivalent width (right panels) for 
the sources in our sample. Where multiple observations exist the weighted mean value has been 
adopted. The upper panels show the values obtained when the line is assumed to be narrow relative 
to the instrumental response (Table ^), whereas the lower plots show the values when the line width 
is left as a free parameter (Table |3|). In both cases, the mean line energy is found to be very close 
to the 6.4 keV expected of neutral iron, even being consistent with a mild redshift. The equivalent 
widths in the broad-line fits are typically ~ 50 per cent larger than for a narrow line. 
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Fig. 2. — Histogram of the line widths obtained from the broad-gaussian fit. In the cases where 
no significant improvement is obtained when a^a is a free parameter, the point has been added 
to the lowest bin. In 14 of the 18 sources the line is resolved, with < (TKa >= 0.43 ± 0.12, which 
corresponds to a FWHM~ 50000 km s -1 . This implies highly relativistic motions if the width is 
due to velocity broadening. 
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Fig. 3. — The profiles of the Fe Ka emission lines in our sample sources, measured by the SIS only. 
The continuum has been estimated by fitting a power-law (with Galactic absorption) to the 3-5 
and 7-10 keV SIS data, and rebinning the resultant residuals. Most show significant broadening 
beyond the instrumental resolution and some have the characteristic profile expected if the line 
is produced close to a central black hole: an asymmetry skewed to the red. All plots have been 
rescaled to a fixed fraction of the continuum at 6.4 keV. That energy (in the rest frame) is marked 
by a vertical dotted line in each plot. The first plot (Mrk 335) shows a horizontal scale in velocity 
space (corresponding to 40,000 km s _1 and a vertical bar which indicates a 10 per cent variation 
from the continuum. 
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Fig. 4. — Mean line profiles for the sample. These were produced by transforming the data/model 
ratios for each sources into the rest-frame of each sources, rebinning the resultant residuals, and 
converting to flux space assuming the mean of the best-fit continua. a) shows the results for the 
whole sample, while b) shows the profile with the two previously-reported skewed lines, MCG-6- 
30-15 and NGC 4151, removed. The two plots are clearly remarkably similar. The solid line in a) 
is a double-gaussian fit to the profile. The narrower component peaks at 6.4 keV with a width of 
a ~ 0.1 keV. The broader component is strongly redshifted, with a centroid energy of 6.1 keV and 
a width of 0.7 keV. It also carries the bulk (~ 75 per cent) of the flux. The dotted line in b) is the 
spectral response of the ASCA SIS detectors at 6.4 keV. 
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Fig. 5. — Data/model residuals for a) our ASCA Seyfert 1 sample (top panel, used to create Fig. |]a); 
the Coma Cluster (middle panel) . The continuum here has been fitted with a bremsstrahlung model, 
excluding the 5-7 keV region; and the Dark Earth. Here the continuum has been modeled with a 
power-law excluding regions contaminated by lines, the most prominent of which are iron Ka and 
nickel Ka. Although all three panels show very clear evidence for emission lines, only the Seyfert 
1 ensemble shows a red wing. 



-45 - 



00 



CO 



CD 
£5 



Pi 



c\2 - 




20 40 60 80 
Disk inclination (i) 



00 



CD 



S-H 

CD 
XI 

3 



w - 



<q>=2.5±0.4 



2 4 6 

Emissivity index (q) 



Fig. 6. — Histograms of the parameters from the disk line fits to the ASCA data (see Table £|). a) 
inclination i, b) index of the emissivity function q (see text). The spectra show a significant spread 
of q, which parameterizes the geometry of the system, indicating that this may change from source 
to source. 
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Fig. 7. — Summed model profiles for the sample. The solid line shows the pure disk-line model 
(Table ||) which is very similar in shape to the data shown in Fig |]. The dashed line shows the 
effects of a narrow component. With the resolution of the ASCA detectors, it is difficult to state 
unequivocally whether or not such a component is present. However, including a narrow line has 
little effect on the disk line profile. 
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Fig. 8. — Change in the disk line parameters when a narrow Gaussian, with energy fixed at 
6.4 keV, is introduced into the fit. There is a small reduction in the equivalent width (top panel) 
corresponding roughly to the equivalent width of the narrow line (with a mean ~ 30 eV). i (middle 
panel) is largely unchanged by the narrow component and the mean q (bottom panel) increases 
insignificantly, to 2.8. 
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Fig. 9.— Residual profiles for NGC 3516 (top panel) and NGC 3227 (bottom panel). Data from all 
four instruments have been combined in bins of 0.25 keV width. The best fit Schwarzschild (solid) 
and Kerr (dotted) models are also shown. 
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Fig. 10. — Histograms of a) the photon index of the power law, T3_io, and b) the line equivalent 
width, Wko from disk line-plus-reflection fits to the ASCA data (see Table |6|). The inclination 
and emissivity index show negligible change when reflection is included. In contrast the power-law 
index is steepened by Ar = 0.12 and the mean equivalent width reduces by ~ 50 eV when the 
effects of reflection are accounted for. 



